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A  SENSITIVITY  ANALYSIS  OF 
UPWELLING  RADIANCES  IN  THE  3-5  jun  REGION 
AS  A  FUNCTION  OF  ATMOSPHERIC  CONDITIONS 


1  INTRODUCTION 

There  is  considerable  interest  in  measuring  atmospheric  parameters  from  remote 
sensing  platforms.  For  example,  meteorological  satellites  routinely  measure  atmospheric 
parameters  for  weather  forecasting  purposes  as  well  as  to  study  the  long  term  characteris¬ 
tics  of  the  atmosphere.  Additionally,  remote  sensing  platforms  are  used  to  study  surface 
conditions.  In  these  cases,  the  atmospheric  portion  of  the  signal  must  be  removed  in  order 
to  deduce  the  information  about  the  surface  or  other  features  being  studied. 

Under  the  direction  of  the  Remote  Sensing  Branch  of  the  Naval  Research  Labora¬ 
tory  (NRL),  SPARTA,  Inc.  calculated  upwelling  ground-to-space  radiances  and  kernel 
functions  for  model  atmospheres  under  different  viewing  and  atmospheric  conditions. 
These  sensitivity  calculations  were  then  used  by  NRL  personnel  in  a  feasibility  study  which 
examined  the  possibility  of  obtaining  information  about  atmospheric  constituents  from  a 
remote  sensing  platform  as  it  looks  through  the  atmosphere  at  a  source  of  a  known  black- 
body  temperature. 

Chapter  2  of  this  report  describes  the  general  methodology  for  the  calculations. 
Chapter  3  gives  the  scope  of  the  sensitivity  calculations  for  the  study.  Finally,  Chapter  4 
summarizes  the  results  of  the  sensitivity  study  and  presents  possible  areas  for  future  work. 


2  GENERAL  METHODOLOGY 
2.1  General  Computational  Approach 

The  upwelling  radiance  calculations  were  performed  using  MODTRAN1,  an  at¬ 
mospheric  transmission  and  radiance  code  developed  by  the  Geophysics  Directorate  of  the 
Phillips  Laboratory  (PL/GP).  For  reference  purposes,  a  “default”  scenario  was  first  estab¬ 
lished  in  terms  of  the  MODTRAN  input  parameters.  Table  1  lists  the  parameters  that 
constituted  the  default  scenario.  Note  that  there  is  no  default  surface  blackbody  tempera¬ 
ture.  Instead,  upwelling  radiances  were  always  calculated  for  a  wide  range  of  surface 
blackbody  temperatures.  The  reason  for  this  approach  was  twofold.  First,  the  present 
study  assumed  the  surface  blackbody  temperature  to  be  a  known,  but  adjustable,  parame¬ 
ter.  Second,  the  range  of  surface  blackbody  temperature  causes  the  surface  emission  com¬ 
ponent  to  dominate  the  total  upwelling  radiance  in  the  3  -  5  pm  region  (see  Section  2.3). 

Departures  from  the  default  scenario  were  imposed  in  a  mutually  exclusive  man¬ 
ner.  That  is,  a  particular  MODTRAN  input  parameter  was  varied  while  the  other  parame¬ 
ters  were  fixed  at  their  default  values.  Although  MODTRAN  permits  multiple  runs  within 
a  given  input  file,  the  decision  was  made  to  use  a  separate  input  file  for  each  perturbed 
scenario.  This  approach  made  it  easier  to  create  the  input  files  as  well  as  post  process  the 
output  files. 

In  each  MODTRAN  input  file,  the  atmospheric  constituent  profiles  were  simulated 
with  the  user-defined  atmosphere  option  (MODEL=7  on  Card  1)  so  that  water  vapor  and 
carbon  dioxide  amounts  could  be  perturbed.  Additionally,  the  IEMSCT  parameter  on 
Card  1  was  set  to  2  so  that  radiances  as  a  function  of  wavenumber  were  included  in 
MODTRANs  TAPE7  output.  The  TAPE 7  output  contains 

1.  Frequency  (cm-1) 

2.  Total  transmittance  (-) 

3.  Atmospheric  radiance  (Watts  cm^sr^cm*1)'1) 

4.  Path  scattered  radiance  (Watts  cm'2sr'1(cm'1)'1) 

5.  Single  scattered  radiance  (Watts  cm’?sr1(cm‘1)“1) 

6.  Total  ground  reflected  radiance  (Watts  cm'2sr1(cm'1)'1) 

7.  Direct  reflected  radiance  (Watts  cm'2sr* ^cm*1)*1) 

8.  Total  radiance  (Watts  cm^sr’Hcm*1)*1) 

and  other  parameters  not  used  in  this  study.  For  each  MODTRAN  run,  the  calculations 
were  performed  over  the  spectral  region  1,920  to  3,330  cm*1  (3.0  -  5.21  |im),  which  rep¬ 
resents  the  full  bandpass  for  a  hypothetical  sensor.  The  spectral  radiances  were  printed  out 
in  steps  of  1  cm*1  with  a  resolution  given  by  a  1  cm*1  rectangular  slit  function,  full-width- 
at-half-maximum  (FWHM).  The  spectral  radiances  in  TAPE7  were  then  post  processed  to 
obtain  upwelling  radiances  for  the  full  sensor  bandpass  plus  32  finer  sensor  bandpasses. 
Table  2  lists  the  bandpasses  considered  in  this  study. 


1  Berk,  A.,  Bernstein,  L.S.,  and  Robertson,  D.C.  (1989)  “MODTRAN:  A  Moderate 
Resolution  Model  for  LOWTRAN7,”  Air  Force  Geophysics  Laboratory,  Hanscom  AFB, 
MA,  AFGL-TR-89-0 1 22,  ADA  214337. 
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Table  1.  MODTRAN  Input  Parameters  Used  to  Describe  the  Default  Scenario.  Radiance 
calculations  were  performed  for  a  range  of  surface  blackbody  temperatures  including  1  K 
(atmospheric  radiance  only),  280  K,  290  K,  330  K,  380  K,  450  K,  and  480  K 


INPUT 

PARAMETER 

DEFAULT 

VALUE 

Model  Atmosphere 

Midlatitude  Summer 

Water  Vapor 

Midlatitude  Summer  Profile 

Carbon  Dioxide 

360ppmv 

Surface  Blackbody  Temperature 

None  (see  figure  caption) 

Boundary  Layer  Aerosol 

Rural 

Surface  Visibility 

23  km 

Cloud  Condition 

Clear 

Sensor  Altitude 

100  km 

Surface  Altitude 

0km 

Surface  Albedo 

0.0  (i.e. ,  blackbody) 

Zenith  Viewing  Angle 

Nadir 

Azimuth  Viewing  Angle 

Looking  North 

Julian  Day  of  Year 

180 

Time  of  Day 

100  Local  Time 

(i.e.,  night,  no  solar  contribution) 

Spectral  Region 

1,920  -  3,330  cm'1 
(see  Table  2  for  bandpasses) 

2.2  Representative  Ground-to-Space  Transmissions 

Figure  1  shows  the  ground-to-space  spectral  transmission  for  the  default  scenario 
(see  Table  1).  Figure  1  also  plots  the  general  locations  of  important  molecular  absorbers. 
To  show  overall  features,  the  transmissions  in  Figure  1  are  degraded  to  20  cm*1  resolu¬ 
tion.  Figure  2  shows  the  water  vapor  transmission  curve  (at  1  cm'1  resolution)  along  with 
the  radiance  curves  for  280  K  and  450  K  blackbodies.  Figure  2  indicates  that  since  the 
blackbody  curves  are  spectrally  shifted  according  to  Wien’s  Law,  the  surface  blackbody 
temperature  can  impact  the  relative  radiance  contributions  in  the  sensor  bandpasses.  Fig¬ 
ure  2  also  demonstrates  that  different  surface  blackbody  temperatures  can  be  used  to  focus 
on  different  regions  of  water  vapor  absorption. 
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Table  2.  Sensor  Bandpasses  Considered  in  This  Study 


BAND 

WAVELENGTH 

WAVENUMBER 

NUMBER 

him 

(cm**) 

COMMENT 

1 

3.00  -  5.21 

1,920-3,330 

Full  Bandpass 

2 

3.35  -  5.00 

2,000  -  2,985 

Truncated  Full  Bandpass 

3 

3.40  -  5.00 

2,000  -  2,941 

h 

4 

3.45  -  5.00 

2,000  -  2,898 

tt 

5 

3.50  -  5.00 

2,000  -  2,857 

ii 

6 

3.35  -  5.05 

1,980-2,985 

Truncated  Full  Bandpass 

7 

3.40  -  5.05 

1,980-2,941 

ii 

8 

3.45  -  5.05 

1,980-2,898 

it 

9 

3.50  -  5.05 

1,980-2,857 

ii 

10 

3.35-5.10 

1,960-2,985 

Truncated  Full  Bandpass 

11 

3.40-5.10 

1,960-2,941 

ii 

12 

3.45-  5.10 

1,960-2.898 

ii 

13 

3.50-5.10 

1,960-2,857 

ii 

14 

3.35-5.15 

1,941  -2,985 

Truncated  Full  Bandpass 

15 

3.40-5.15 

1,941  -  2,941 

it 

16 

3.45-5.15 

1,941  -2,898 

•i 

17 

3.50-5.15 

1,941  -2,857 

it 

18 

3.35  -  4.00 

2,500  -  2,985 

Low  Bandpass 

19 

3.40  -  4.00 

2,500  -  2,941 

ii 

20 

3.45  -  4.00 

2,500  -  2,898 

ii 

21 

3.50  -  4.00 

2,500  -  2,857 

ii 

22 

3.35  -  4.05 

2,469  -  2,985 

Low  Bandpass 

23 

3.40  -  4.05 

2,469  -  2,941 

ii 

24 

3.45  -  4.05 

2,469  -  2,898 

n 

25 

3.50  -  4.05 

2,469  -  2,857 

ii 

26 

3.35-4.10 

2,439  -  2,985 

Low  Bandpass 

27 

3.40-4.10 

2,439  -  2,941 

ii 

28 

3.45-4.10 

2,439  -  2,898 

ii 

29 

3.50-4.10 

2,439  -  2,857 

ii 

4.50  -  5.00 

2,000  -  2,222 

High  Bandpass 

4.50  -  5.05 

1,980-2,222 

ii 

4.50-5.10 

1,960-2,222 

ii 

4.50-5.15 

1,941  -2,222 

it 

'  ■m* 
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3.0  3.2  3.4  3.6  3.8  4.0  4.2  4.4  4.6  4.8  5.0  5.2 

Wavelength  (Mm) 

Figure  1.  Ground-to-Space  Transmission  Versus  Wavelength  for  the  Default  Scenario.  In 
order  to  show  overall  features,  the  transmissions  are  degraded  to  20  cm'1  resolution 


Wavelength  (jum) 

Figure  2.  Blackbody  Spectral  Radiance  and  Ground-to-Space  Water  Vapor  Transmission 
Versus  Wavelength.  The  calculations  represent  the  default  conditions  in  Table  1.  The 
transmission  curve  was  generated  from  a  10  cm'1  grid  of  values  at  1  cm'1  resolution 
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To  help  establish  the  impact  of  water  vapor  on  the  upwelling  radiance,  transmis¬ 
sion  calculations  similar  to  those  in  Figure  1  were  performed  using  extreme  water  vapor 
profiles.  The  extreme  water  vapor  profiles  (and  the  accompanying  profiles  of  pressure, 
temperature,  etc.)  were  represented  by  the  Tropical  and  Subarctic  Winter  AFGL  model 
atmospheres.2  Figure  3  shows  the  ground- to-space  water  vapor  transmission  for  Tropical 
and  Subarctic  Winter  model  atmospheres,  plus  their  ratio.  To  show  overall  features,  the 
transmissions  are  degraded  to  20  cm-1  resolution.  Figure  3  indicates  that  the  greatest 
changes  in  water  vapor  transmission  are  between  3.0  and  3.4  pm  and  between  4.5  and 
5.21  pm. 


Wavelength  (Mm) 

Figure  3.  Ground-to-Space  Water  Vapor  Transmission  Versus  Wavelength  for  Tropical 
and  Subarctic  Winter  Model  Atmospheres.  The  ratio  of  the  water  vapor  transmittances  is 
also  shown.  In  order  to  show  overall  features,  the  transmissions  are  degraded  to  20  cm*1 
resolution 


2.3  Representative  Impact  of  the  Surface  Blackbody  Temperature 

Since  this  study  involved  a  wide  range  of  surface  blackbody  temperatures,  a  set  of 
preliminary  calculations  established  that  the  upwelling  radiance  depends  strongly  on  the 
surface  blackbody  temperature.  For  example,  Figure  4  shows  the  upwelling  radiance  for 
surface  blackbody  temperatures  of  1,  280,  330,  380,  and  450  K.  The  default  scenario  is 


2  Anderson,  G.P.,  Clough,  S.A.,  Kneizys,  F.X.,  Chetwynd,  J.H.,  and  Shettle,  E.P.  (1986) 
“AFGL  Atmospheric  Constituent  Profiles  (0-120  km),”  Air  Force  Geophysics  Labora¬ 
tory,  Hanscom  AFB,  MA  01731,  AFGL-TR-86-01 10, 15  May,  ADA  175173. 
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used  for  these  calculations  and,  to  show  overall  features,  die  radiances  are  degraded  to 
20  cm'1  resolution. 

The  impact  of  the  surface  blackbody  temperature  on  different  sensor  bandpasses 
was  also  examined.  Figure  5  shows  integrated  radiances  for  33  sensor  bandpasses  (see 
Table  2)  and  different  surface  blackbody  temperatures. 


Wavelength  (Mm) 

Figure  4.  Ground-to-Space  Spectral  Radiance  Versus  Wavelength  For  Different  Surface 
Blackbody  Temperatures.  Radiances  for  a  1  K  blackbody  represent  atmospheric  contribu¬ 
tions  only.  The  calculations  represent  the  default  scenario  in  Table  1.  In  order  to  show 
overall  features,  radiances  are  degraded  to  20  cm*1  resolution 
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Figure  5.  Integrated  Ground-to-Space  Radiances  for  33  Sensor  Bandpasses  and  Different 
Surface  Blackbody  Temperatures.  Integrated  radiances  for  a  1  K  blackbody  represent  at¬ 
mospheric  contributions  only.  The  calculations  represent  the  default  scenario  in  Table  1 
and  the  sensor  bandpasses  in  Table  2 


3  SCOPE  OF  THE  SENSITIVITY  ANALYSES 

A  major  goal  of  this  study  was  to  examine  how  variations  in  the  viewing  and  at¬ 
mospheric  conditions  impact  the  upwelling  radiance.  To  do  this,  upwelling  radiances  were 
calculated  for  “typical”  environmental  conditions  and  different  viewing  geometries.  For 
reference.  Table  3  gives  the  values  that  were  used  for  the  viewing  and  atmospheric  pa¬ 
rameters. 

In  the  current  study,  an  independent  sensitivity  analysis  was  performed  for  each 
viewing  and  atmospheric  parameter.  That  is,  while  one  parameter  was  varied,  the  other 
parameters  were  fixed  at  their  default  values  (unless  stated  otherwise).  In  each  analysis, 
upwelling  radiances  were  calculated  for  the  full  range  of  surface  blackbody  temperatures 
and  for  each  sensor  bandpass  listed  in  Table  2.  Information  about  each  analysis  is  given  in 
the  following  sections.  Sample  results  are  shown  in  Chapter  4.  Because  the  sensitivity 
analyses  created  large  amounts  of  data,  only  a  portion  of  the  results  are  presented.  The  full 
set  of  tabulated  radiances  can  be  obtained  from  SPARTA,  Inc. 
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Table  3.  Summary  of  the  Sensitivity  Analyses  in  This  Study.  The  table  also  lists  the  per¬ 
turbed  values  and,  where  possible,  the  default  values  are  denoted  by  an  asterisk 


PARAMETER 

UNITS 

VALUES 

Water  Vapor  Mixing  Ratio 

ppmv 

Thirteen  profiles  (see  Section  3.1) 

Carbon  Dioxide  Mixing  Ratio 

ppmv 

330.0,  359.9, 360.0*,  360.1, 390.0 

Sensor  Altitude 

km 

100*,  25, 15 

Zenith  Viewing  Angle  ' 

deg 

Nadir*  and  30, 45,  60  Off  Nadir 

Surface  Visibility 

km 

50,  23*,  5 

Cloud  Condition 

- 

Clear*,  Thin  Cirrust 

Solar  Contribution 

Local  Time 

0100  (Nighttime)*,  1200  hrs 

f  The  range  of  extinction  coefficients  for  the  thin  cirrus  was  taken  as  0.014, 0.028,  and 
0.056  km'l  at  0.55  pm 


3.1  Variations  in  Water  Vapor  Profile 

A  major  requirement  of  this  study  was  to  see  how  different  water  vapor  profiles 
impact  the  upwelling  radiance  since  it  is  a  highly  variable  atmospheric  quantity.  Therefore, 
the  water  vapor  amount  was  perturbed  extensively  to  satisfy  this  requirement 

To  easily  represent  as  many  realistic  water  vapor  profiles  as  possible,  perturbed 
water  vapor  amounts  were  described  as  “percent  perturbations”  from  the  reference  profile. 
As  an  example,  suppose  the  reference  water  vapor  mixing  ratio  amounts  were  20,000 
ppmv  and  10,000  ppmv  at  altitudes  of  0  km  and  2  km,  respectively.  For  a  -20%  perturba¬ 
tion,  the  perturbed  water  vapor  mixing  ratio  amount  equals 

20,000  -  20,000  x  (.20)  =  16,000  ppmv  at  0  km 
and 

10,000  -  10,000  x  (.20)  =  8,000  ppmv  at  2  km. 

In  this  study,  twelve  perturbed  water  vapor  mixing  ratio  profiles  were  constructed 
from  a  reference  water  vapor  profile.  The  percent  perturbation  amounts  for  water  vapor 
are  listed  in  Table  4.  For  the  reference  water  vapor  profiles.  Midlatitude  Summer  and 
Winter  conditions  were  used  along  the  accompanying  profiles  of  pressure,  temperature, 
etc.  Also,  the  sensitivity  analysis  considered  surface  blackbody  temperatures  of  1,  280, 
330,  380,  and  450K,  and  zenith  viewing  angles  of  nadir,  and  30,  45,  and  60  degrees  off 
nadir.  Thus,  this  investigation  consisted  of  13x2x5x4=520  different  scenarios.  Other  input 
parameters  were  fixed  at  their  default  values. 
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Table  4.  Summary  of  the  Sensitivity  Analysis  for  Water  Vapor.  The  perturbations  are 
performed  by  applying  a  percent  variation,  relative  to  the  reference  amount  in  each  layer, 
over  the  listed  altitude  region 


PERTURBATION 

NUMBER 

PERCENT 

PERTURBATION 

PERTURBATION 
ALTITUDE  REGION 

1 

-20% 

0-  100  km 

2 

-10% 

0-  100  km 

3 

-10% 

0  -  5  km  (0%  elsewhere) 

4 

-10% 

6  -  15  km  (0%  elsewhere) 

5 

-5% 

6-15  km  (-1%  elsewhere) 

6 

-1% 

0-  100  km 

7 

0%  (Reference) 

0-  100  km 

8 

1% 

0-  100  km 

9 

5% 

6-15  km  (1%  elsewhere) 

10 

10% 

6  - 15  km  (0%  elsewhere) 

11 

10% 

0  -  5  km  (0%  elsewhere) 

12 

10% 

0-  100  km 

13 

20% 

0-  100  km 

3.2  Variations  in  Carbon  Dioxide  Profile 

The  impact  of  the  carbon  dioxide  amount  on  the  upwelling  radiance  was  also  in¬ 
vestigated.  This  analysis  was  done  because  the  bandpass  regions  contain  a  very  strong 
C02  absorption  band  between  about  4. 1  and  4.5  (im  and  because  the  amount  of  C02  in  the 
atmosphere  varies  seasonally.  Thus,  it  was  surmised  that  the  radiance  contributions  from 
the  edges  of  the  band  might  be  affected. 

Using  a  default  value  of  360  ppmv  throughout  the  atmosphere,  C02  amounts  were 
varied  with  ±10%  and  ±1%  perturbations.  The  model  atmosphere  profiles  for  Midlatitude 
Summer  and  Winter  conditions  were  used  for  the  remaining  atmospheric  constituents.  The 
sensitivity  analysis  considered  surface  blackbody  temperatures  of  290,  330,  and  480  K. 
Thus,  this  investigation  consisted  of  5x2x3=30  different  scenarios.  Other  input  parameters 
were  fixed  at  their  default  values. 


3.3  Variations  in  Sensor  Altitude 

The  impact  of  the  sensor  altitude  on  the  upwelling  radiance  was  investigated.  This 
analysis  was  performed  to  see  if  the  presence  of  middle  and  upper  atmosphere  constitu¬ 
ents,  such  as  ozone  and  stratospheric  aerosols,  affect  the  upwelling  radiance. 

Upwelling  radiances  were  calculated  for  sensor  altitudes  of  15,  25,  and  100  km. 
The  sensitivity  analysis  considered  surface  blackbody  temperatures  of  290, 330,  and  480  K 
and  model  atmosphere  profiles  for  Midlatitude  Summer  and  Winter  conditions.  Thus,  this 
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investigation  consisted  of  3x3x2=18  different  scenarios.  Other  input  parameters  were 
fixed  at  their  default  values. 


3.4  Variations  In  Zenith  Viewing  Angle 

The  impact  of  the  zenith  viewing  angle  on  the  upwelling  radiance  was  investigated. 
This  analysis  was  performed  because  the  sensor  may  be  required  to  measure  the  upwelling 
radiance  for  off-nadir  viewing  conditions.  Since  the  path  length  through  the  atmosphere 
depends  on  the  zenith  viewing  angle,  it  was  surmised  that  the  upwelling  radiance  would  be 
affected. 

Upwelling  radiances  were  calculated  for  nadir  viewing  conditions  and  for  30,  45, 
and  60  degrees  off  nadir.  The  sensitivity  analysis  considered  surface  blackbody  tempera¬ 
tures  of  290,  330,  and  480  K  and  model  atmosphere  profiles  for  Midlatitude  Summer  and 
Winter  conditions.  Thus,  this  investigation  consisted  of  4x3x2=24  different  scenarios. 
Other  input  parameters  were  fixed  at  their  default  values. 


3.5  Variations  in  Surface  Visibility 

The  impact  of  the  surface  visibility  on  the  upwelling  radiance  was  investigated. 
This  analysis  was  performed  because  the  boundary  layer  aerosols  attenuate  the  emitted 
radiance  from  the  surface  blackbody  source  plus  they  scatter  and  absorb  the  atmospheric 
path  radiance. 

Upwelling  radiances  were  calculated  for  surface  visibilities  of  5,  23,  and  50  km. 
The  sensitivity  analysis  considered  surface  blackbody  temperatures  of  290,  330,  and  480K 
and  model  atmosphere  profiles  for  Midlatitude  Summer  and  Winter  conditions.  Thus,  this 
investigation  consisted  of  3x3x2=18  different  scenarios.  Other  input  parameters  were 
fixed  at  their  default  values. 


3.6  Variations  in  Cloud  Conditions 

The  impact  of  the  cloud  cover  on  the  upwelling  radiance  was  investigated.  This 
analysis  was  performed  because  thin  cirrus  might  attenuate  the  emitted  radiance  from  the 
surface  blackbody  source  plus  they  also  might  contribute  to  the  atmospheric  radiance. 
Low  and  middle  level  clouds,  such  as  cumulus  and  altocumulus,  were  not  considered  be¬ 
cause  they  are  optically  thick  and  completely  block  out  the  surface  emission  term. 

Upwelling  radiances  were  calculated  for  clear  skies  and  for  three  thin  cirrus  types. 
To  account  for  their  highly  variable  nature,  the  three  thin  cirrus  were  assigned  extinction 
coefficients  (at  0.55  |im)  equal  to  0.014,  0.028,  and  0.056  km*1.  The  cirrus  thickness  was 
fixed  at  0.2  km.  The  sensitivity  analysis  considered  surface  blackbody  temperatures  of  1, 
280,  330,  380,  and  480K  and  model  atmosphere  profiles  for  Midlatitude  Summer  and 
Winter  conditions.  Also,  because  the  cirrus  optical  depth  increases  greatly  when  viewed 
off  nadir,  the  zenith  viewing  angle  was  varied  as  nadir  and  30,  45,  and  60  degrees  off  na¬ 
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dir.  Thus,  this  investigation  consisted  of  4x5x2x4=160  different  scenarios.  Other  input 
parameters  were  fixed  at  their  default  values. 


3.7  Variations  in  Solar  Contribution 

The  impact  of  the  solar  contribution  on  the  upwelling  radiance  was  investigated. 
This  analysis  was  performed  because  a  small  portion  of  the  solar  spectrum  overlaps  the 
sensor  bandpasses  near  3.0  |im.  Although  ground-reflected  solar  radiance  is  zero  (because 
the  sensor  is  looking  at  a  known  blackbody),  this  analysis  assesses  solar  scattering  off 
aerosols  and  molecules  along  the  slant  path. 

To  evaluate  the  solar  contribution,  upwelling  radiances  were  calculated  for  night¬ 
time  and  1200  local  time.  The  sensitivity  analysis  considered  surface  blackbody  tempe 
tures  of  290,  330,  and  480K  and  model  atmosphere  profiles  for  Midlatitude  Summer 
Winter  conditions.  Thus,  this  investigation  consisted  of  2x3x2=12  different  scenarios 
Other  input  parameters  were  fixed  at  their  default  values. 


3.8  Kernel  Calculations 

In  addition  to  upwelling  radiance,  the  monochromatic  and  broadband  kernel  func¬ 
tions,  k(z,v)  and  K(z)  respectively,  were  calculated.  Physically,  the  kernel  function  de¬ 
scribes  the  rate  of  change  of  transmission  as  a  function  of  altitude  and  can  be  used  to  un¬ 
derstand  what  regions  of  the  atmosphere  contribute  the  most  to  the  upwelling  radiance. 


3.8.1  Governing  Equations 


The  monochromatic  kernel  function  is  given  by 


fc(z,v) 


dx(z,v) 

dz 


(1) 


where  z  and  v  are  altitude  and  frequency,  respectively,  and  x(z,v)  is  the  monochromatic 
transmission  along  the  slant  path.  The  broadband  kernel  function  is  given  by 

Pj k(z,v)b(z,v)dv 

K(z)  =  ^~p~ -  (2) 

J  b{z,v)d» 


where  b(z,v)  is  the  monochromatic  Planck  radiance  and  V;  and  V2  are,  respectively,  the 
lower  and  upper  limits  of  the  sensor  bandpass. 

The  monochromatic  and  broadband  kernel  functions  are  relevant  to  this  study  be¬ 
cause  the  atmospheric  contributions  to  the  upwelling  path  radiance  can  be  evaluated  in 
terms  of  them.  In  a  layered  atmosphere  in  which  scattering  is  negligible  compared  to  ther¬ 
mal  emission,  the  monochromatic  path  radiance,  ifz^v),  is  given  by 
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Kz,,v)  =  J*  k(z,v)b(z,v)dz 

(3) 

where  zs  is  the  sensor  altitude.  Note  that  the  expression  for  i(zs,v)  does  not  include  the 
surface  blackbody  component.  Also,  the  expression  for  i(zs>v)  is  approximate  and  does  not 
form  the  basis  for  the  MODTRAN  radiance  calculations  in  this  study.  Similarly,  it  can  be 
shown  that  the  broadband  path  radiance,  I(zs),  is  given  by 

I(zs)  =  fiK(z)B(z)dz 

(4) 

where  B(z)  is  the  band-integrated  Plank  radiance 

B(z)  =  f  *b(z,\)cN 

(5) 

Eqs.  3  and  4  indicate  that  the  profiles  of  the  monochromatic  and  broadband  kernel  func¬ 
tions  can  be  used  to  estimate  the  relative  amounts  of  thermal  emission  from  the  altitude 
regions  along  the  slant  path.  Specifically,  the  maxima  of  k(z,v)  and  K(z)  indicate  the  alti¬ 
tudes  that  contribute  the  most  to  the  upwelling  radiance. 


3.8.2  Method  for  Calculating  Kernel  Functions 

The  method  for  calculating  monochromatic  and  broadband  kernel  functions  re¬ 
quires  extensive  post  processing  of  the  standard  MODTRAN  output.  Thus,  the  mono¬ 
chromatic  and  broadband  kernel  functions  were  only  calculated  for  a  few  representative 
scenarios.  Table  5  lists  the  scenarios  that  were  considered.  To  automate  the  process, 
monochromatic  kernel  functions  were  calculated  using  SENTRAN7.3  The  general  proce¬ 
dure  was  to 

1.  Define  a  scenario  in  SENTRAN7's  Edit  Module.  (This  module  mimics  the  flow 
through  MODTRAN’s  input  deck.) 

2.  Perform  MODTRAN  calculations  for  the  scenario. 

3.  Enter  SENTRAN7's  Graph  and  Analyze  Module  and  select  the  option  that  calcu¬ 
lates  monochromatic  kernel  functions.  (The  Graph  and  Analyze  Module  allows  users 
to  plot  and  manipulate  the  standard  MODTRAN  output  The  option  to  calculate 
monochromatic  kernel  functions  is  a  standard  feature  in  SENTRAN7.) 

4.  Instruct  SENTRAN7  to  write  the  monochromatic  kernel  functions  to  a  data  file. 

5.  Select  the  option  in  the  Graph  and  Analyze  Module  that  calculates  the  product  of 
the  monochromatic  kernel  and  blackbody  functions.  (This  option  is  standard  in 
SENTRAN7.) 


3  Longtin,  D.R.,  DePiero,  N.L.,  Pagliughi,  F.P.,  and  Hummel,  J.R.  (1991)  “SENTRAN7: 
The  Sensitivity  Analysis  Package  for  LOWTRAN7  and  MODTRAN,”  Phillips  Labora¬ 
tory,  Hanscom  AFB,  MA,  PL-TR-91 -2290(D) ,  ADA  251595. 
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6.  Instruct  SENTRAN7  to  write  the  product  of  the  monochromatic  kernel  and  black- 
body  functions  to  a  data  file. 

7.  Calculate  the  broadband  kernel  functions  off-line  using  the  data  file  in  Step  6. 

This  procedure  was  repeated  for  each  scenario  in  Table  5. 


Table  5.  Scenarios  for  the  Monochromatic  and  Broadband  Kernel  Function  Calculations. 
Water  vapor  profiles  are  based  on  the  model  atmospheres.  The  “full  bandpass”  covers 
3-5.21  pm,  the  “truncated  full  bandpass”  covers  3.45  -  5.15  pm,  and  the  “high  bandpass” 
covers  4.5  -  5.15  pm.  Note  that  the  kernel  functions  only  consider  the  atmospheric  path  so 
surface  blackbody  temperatures  are  not  required 


ZENITH 

MODEL 

SCENARIO 

VIEWING 

ATMOSPHERE 

SENSOR 

NUMBER 

ANGLE 

PROFILE 

BANDPASS 

1 

Nadir 

Midlatitude  Summer 

Full  Bandpass 

2 

45°  off  Nadir 

Midlatitude  Summer 

it 

3 

Nadir 

Midlatitude  Winter 

it 

4 

Nadir 

Tropical 

If 

5 

Nadir 

Subarctic  Winter 

H 

6 

Nadir 

Midlatitude  Summer 

Truncated  Full  Bandpass 

7 

45°  off  Nadir 

Midlatitude  Summer 

II 

8 

Nadir 

Midlatitude  Winter 

if 

9 

Nadir 

Tropical 

if 

10 

Nadir 

Subarctic  Winter 

it 

11 

Nadir 

Midlatitude  Summer 

High  Bandpass 

12 

45°  off  Nadir 

Midlatitude  Summer 

ti 

13 

Nadir 

Midlatitude  Winter 

it 

14 

Nadir 

Tropical 

it 

15 

Nadir 

Subarctic  Winter 

if 

4  RESULTS  FROM  THE  SENSITIVITY  ANALYSES 

Tabulated  values  for  all  the  sensitivity  analyses  have  been  assembled,  but  they  are 
too  extensive  to  include  in  this  report.  The  results  can  be  made  available  to  interested  in¬ 
dividuals.  In  this  chapter,  the  primary  findings  from  the  sensitivity  analyses  are  presented. 
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4.1  Water  Vapor  Analysis 


The  impact  of  water  vapor  amount  on  the  upwelling  radiance  depends  on  where  in 
the  atmosphere  the  perturbations  exist,  the  viewing  angle,  the  sensor  bandpass,  and  the 
model  atmosphere.  Figure  6  shows  sample  upwelling  radiances  for  one  set  of  water  vapor 
profiles  as  a  function  of  the  surface  blackbody  temperature.  In  the  example,  the  conditions 
represent  a  Midlatitude  Summer  atmosphere  and  the  perturbations  were  applied  at  all  alti¬ 
tudes  in  the  atmosphere.  In  Figure  7,  the  differences  in  the  radiances  from  the  perturbed 
water  vapor  profiles  (relative  to  the  default  conditions)  are  shown  as  a  function  of  the 
surface  blackbody  temperature.  Figure  8  displays  the  differences  in  the  radiances  when  the 
perturbations  occur  in  different  regions  of  the  atmosphere. 

A  comparison  of  Figures  7  and  8  suggests  that  water  vapor  perturbations  in  the 
lower  atmosphere  contribute  the  most  to  the  change  in  the  upwelling  radiance.  In  this 
case,  the  lower  atmosphere  is  more  important  because  water  vapor  amounts  are  greater. 
Also,  the  impact  on  the  upwelling  radiance  becomes  greater  when  the  temperature  of  the 
surface  blackbody  source  is  increased.  This  behavior  was  found  to  be  true  for  all  of  the 
model  atmospheres. 


275  300  325  350  375  400  425  450 

Blackbody  Temperature  (K) 

Figure  6.  Upwelling  Radiance  Versus  Surface  Blackbody  Temperature  for  Perturbed 
Profiles  of  Water  Vapor  Mixing  Ratio 
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275  300  325  350  375  400  425  450 

Blackbody  Temperature  (K) 


Figure  7.  Differences  in  Upweiling  Radiance  Versus  Surface  Blackbody  Temperature  for 
Perturbed  Profiles  of  Water  Vapor  Mixing  Ratio 


Blackbody  Temperature  (K) 

Figure  8.  Differences  in  Upweiling  Radiance  Versus  Surface  Blackbody  Temperature  for 
Perturbed  Water  Vapor  Profiles  in  Different  Regions  of  the  Atmosphere.  In  this  example, 
the  perturbations  were  applied  in  0  -  5  km  region  (LL)  and  in  the  6  - 15  km  region  (ML) 
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4.2  Other  Sensitivity  Analyses 

Figure  9  shows  representative  upwelling  radiances  as  a  function  of  surface  black- 
body  temperature  for  different  viewing  angles.  Figure  10  shows  representative  upwelling 
radiances  as  a  function  of  surface  blackbody  temperature  for  different  values  of  cirrus  ex¬ 
tinction  coefficient  The  corresponding  differences  from  reference  values  (i.e.,  clear  skies) 
are  shown  in  Figure  1 1 . 

In  order  to  evaluate  the  relative  impact  of  each  perturbed  parameter,  the  percent 
difference  was  calculated  relative  to  the  reference  upwelling  radiance.  Specifically,  the 
percent  difference  was  computed  as 

%  dif  100  *  ^ax™um  Radiance-  Minimum  Radiance) 

Reference  Radiance 

for  each  combination  of  model  atmosphere,  surface  blackbody  temperature,  and  sensor 
bandpass.  The  percent  differences  are  summarized  in  Table  6. 


275  300  325  350  375  400  425  450 

Blackbody  Temperature  (K) 


Figure  9.  Upwelling  Radiance  Versus  Surface  Blackbody  Temperature  for  Different  Ze¬ 
nith  Viewing  Angles 
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Blackbody  Temperature  (K) 

Figure  10.  Upwelling  Radiance  Versus  Surface  Blackbody  Temperature  for  Different  Cir¬ 
rus  Cloud  Extinction 


275  300  325  350  375  400  425  450 

Blackbody  Temperature  (K) 

Figure  11.  Differences  in  Upwelling  Radiance  Versus  Surface  Blackbody  Temperature  for 
Different  Cirrus  Cloud  Extinction 
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Tabic  6.  Percent  Differences  in  Upwelling  Radiance  for  the  Different  Sensitivity  Analyses. 
The  bandpass  names  correspond  to  those  given  in  Table  2 


PERTURBATION 

ANALYSIS 

MODEL 

ATMOSPHERE 

SURFACE 

BLACKBODY 

TEMPERATURE 

(K) 

SENSOR 

BANDPASS 

PERCENT 

DIFFERENCE 

Carbon  Dioxide 

Midlatitude  FW* **  and  SS 

290,  330, 480 

All  Bands 

<1% 

Sensor  Altitude 

Midlatitude  FW  and  SS 

290,  330, 480 

Full  Bands 

1-3% 

Midlatitude  FW  and  SS 

290,  330,  480 

Low  Bands 

<1% 

Midlatitude  FW  and  SS 

290,  330,  480 

High  Bands 

3-5% 

Zenith  Viewing  Angle 

Midlatitude  SS 

290 

Full,  High  Bands 

14-16% 

Midlatitude  SS 

290 

Low  Bands 

5-7% 

Midktitude  FW 

290 

Full  High  Bands 

19-21% 

Midlatitude  FW 

290 

Low  Bands 

9-11% 

Midlatitude  FW  and  SS 

330 

Full  Bands 

24-28% 

Midlatitude  FW  and  SS 

330 

Low  Bands 

12-16% 

Midlatitude  FW  and  SS 

330 

High  Bands 

29-33% 

Midladtude  SS 

480 

Full  Bands 

30-33% 

Midlatitude  SS 

480 

Low  Bands 

19-21% 

Midlatitude  SS 

480 

High  Bands 

48-50% 

Midlatitude  FW 

480 

Full  Bands 

24-27% 

Midlatitude  FW 

480  • 

Low  Bands 

12-15% 

Midlatitude  FW 

480 

High  Bands 

32-36% 

Surface  Visibility 

Midlatitude  SS 

290 

All  Bands 

<1% 

Midlatitude  FW 

290 

All  Bands 

3-6% 

Midladtude  SS 

330 

Full,  High  Bands 

4-6% 

Midlatitude  SS 

330 

Low  Bands 

7-8% 

Midladtude  FW 

330 

All  Bands 

7-9% 

Midladtude  FW  and  SS 

480 

All  Bands 

8-10% 

Cloud  Condition 

Midladtude  SS 

1” 

Full  High  Bands 

<l%t 

Midlatitude  SS 

l" 

Low  Bands 

5-6%  t 

Midladtude  FW 

l” 

Full,  High  Bands 

1-3%  t 

Midladtude  FW 

l” 

Low  Bands 

8-1  l%t 

Midlatitude  FW  and  SS 

280 

All  Bands 

<l%t 

Midlatitude  FW  and  SS 

330,  380, 450 

All  Bands 

1-3%  t 

Solar  Contribution 

Midlatitude  FW  and  SS 

290 

Low  Bands 

1-2% 

Midlatitude  FW  and  SS 

290 

Full,  High  Bands 

<1%  ! 

Midlatitude  FW  and  SS 

330, 480 

All  Bands 

<  i% 

*  FW  =  Fall-Winter,  SS  =  Spring-Summer 

** Atmospheric  contributions  only 
t  All  viewing  angles  considered 
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From  Table  6,  it  can  be  concluded  that  for  the  sensor  bandpasses  investigated  in 
this  study, 

•  The  impact  of  the  temperature  of  the  surface  blackbody  source  is  significant 

•  The  impact  of  the  C02  profile  is  small 

•  The  impact  of  the  sensor  altitude  is  small 

•  The  choice  of  viewing  angle  is  significant  and  a  function  of  the  sensor  bandpass 

•  The  impact  of  surface  visibility  is  small 

•  The  impact  of  solar  contributions  (i.e.  day/night  differences)  is  insignificant 


4.3  Broadband  Kernel  Functions 

Figure  12  shows  representative  broadband  kernel  functions  for  this  study.  In  the 
example  shown,  the  kernel  functions  were  calculated  for  the  full  bandpass  (1,920  - 
3,330  cm ')  and  for  two  model  atmospheres  and  viewing  conditions.  The  kernel  function 
peaks  in  the  lower  atmosphere  in  the  2  -  4  km  range*  indicating  the  greatest  sensitivity  to 
changes  in  atmospheric  transmission.  The  basic  shape  of  the  kernel  function  profile  is,  in 
general,  insensitive  to  the  choice  of  model  atmosphere,  viewing  angle,  and  sensor  band¬ 
pass,  as  shown  in  Figure  13.  However,  the  magnitude  of  the  kernel  function  at  a  given 
altitude  depends  on  the  choice  of  viewing  angle  and  sensor  bandpass. 
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Figure  12.  Sample  Broadband  Kernel  Functions  Versus  Altitude  for  Different  Model  At¬ 
mospheres  and  Viewing  Conditions 
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Figure  13.  Broadband  Kerne'  Functions  Versus  Altitude  for  Different  Sensor  Bandpasses 
and  a  Viewing  Angle  of  (a.)  180°  and  (b.)  135° 
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4.4  Future  Work 


The  concept  of  obtaining  information  about  atmospheric  constituents  through  the 
use  of  elevated  background  temperatures  is  viable  and  research  on  it  should  continue.  For 
water  vapor  applications,  it  is  quite  possible  that  the  upwelling  radiance  in  other  wave¬ 
length  regions  (i.e.,  sensor  bandpasses)  can  be  used  in  conjunction  with  that  in  3-5  pm  re¬ 
gion  as  a  way  to  improve  the  effectiveness  of  the  approach.  These  additional  wavelength 
regions  should  be  located.  Candidate  regions  must  satisfy  a  few  simple  criteria.  For  exam¬ 
ple,  the  transmission  in  the  region  must  be  highly  sensitive  to  changes  in  water  vapor 
amount.  Also,  the  region  must  not  be  opaque  (including  any  effects  from  aerosol  and  mo¬ 
lecular  scattering).  Finally,  the  region  should  be  insensitive  to  changes  in  other  atmos¬ 
pheric  constituents  and  ideally,  the  transmission  for  other  atmospheric  constituents  should 
be  near  1.0. 

The  candidate  wavelength  regions  can  be  found  using  MODTRAN  calculations  for 
Tropical  and  Subarctic  Winter  conditions.  Beside  checking  for  opacity,  maximum  values 
of  a  simple  parameterization  such  as 

P(\)  -  IAThv(X,)/ATot7/£r(^«)I  •  (7) 

can  be  used  to  rapidly  identify  candidate  regions.  In  Eq.  7,  Axm  is  the  change  in  water  va¬ 
por  transmission  for  Tropical  and  Subarctic  Winter  conditions  and  Axomm  is  the  change  in 
transmission  due  to  other  processes  (scattering  and  other  molecular  absorption)  for  Tropi¬ 
cal  and  Subarctic  Winter  conditions.  Next,  the  candidate  regions  should  be  refined  further 
by  calculating  the  corresponding  monochromatic  and  broadband  kernel  functions.  For  re¬ 
mote  sensing  applications  that  require  information  about  water  vapor  in  a  particular  alti¬ 
tude  region,  the  kernel  function  should  reach  its  maximum  value  in  that  region. 

The  current  concept  can  also  be  extended  to  other  types  of  gases  and  molecules  of 
interest.  Possible  applications  would  be  to  locate  or  track  the  presence  of  toxic  fumes  or 
chemical  biological  warfare  (CBW)  agents  near  forest  fires  or  structural  fires  (i.e.,  a  hot 
background).  To  obtain  the  best  wavelength  region  for  the  application,  the  analysis  de¬ 
scribed  above  should  be  performed  for  the  gas  or  molecule  of  interest. 
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